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ABSTRACT 

Mesoscale  vortices  generated  by  western  boundary  currents  are  well 
observed  and  documented,  particularly  in  the  case  of  the  Gulf  Stream 
System.  The  movement  of  these  rings  in  the  region  of  the  Gulf  Stream 
is  well  studied  and  has  been  ascribed  to  the  following  physical 
mechanisms:  (1)  the  beta  effect  on  an  isolated  ring,  (2)  advection  of 
a  ring  in  a  recirculation  regime,  (3)  downstream  advection  of  a  ring 
in  contact  with  a  jet,  and  (4)  vorticity  advection  associated  with  the 
jet  and  eddy  interaction. 

Utilizing  a  two  layer,  nonlinear  primitive  equation  model,  an 
examination  of  eddy  movement  is  conducted,  with  focus  on  eddy/ jet 
interaction.  A  series  of  numerical  experiments  is  performed  in  which 
the  initial  separation  distance  between  eddy  and  jet  is  varied.  The 
model  demonstrates  that  vortex  movement  is  strongly  related  to  the 
proximity  of  the  vortex  to  the  jet.  It  also  is  demonstrated  that 
observed  movement  is  not  solely  dependent  on  the  beta  effect  nor  on 
advection  due  to  recirculation. 
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I.   INTRODUCTION 

A.   THE  GULF  STREAM  SYSTEM 

The  Gulf  Stream  System  is  an  extensively  studied  current  system 
due  both  to  its  unique  dynamical  properties  as  well  as  to  the  wealth 
of  accumulated  data.  It  is  the  strongest  current  system  of  the  North 
Atlantic  subtropical  gyre.  Iselin  (1936)  proposed  terminology  for  a 
Gulf  Stream  System  subdivided  into  three  distinctly  different  currents 
that  is  commonly  accepted  and  will  be  used  in  this  study.  The  first 
current,  the  Florida  Current,  is  the  portion  of  the  Gulf  Stream 
System  that  originates  in  the  Gulf  of  Mexico  with  the  waters  that 
flow  through  the  Yucatan  Channel.  It  passes  through  the  Florida 
Straits  and  travels  northward  along  the  eastern  United  States  to  the 
vicinity  of  Cape  Hatteras.  At  this  juncture  it  turns  eastward  and 
leaves  the  continental  slope  region  towards  deeper  waters.  Here  it  is 
known  as  the  Gulf  Stream  as  it  flows  towards  the  Grand  Banks  region. 
It  is  within  this  region  that  the  stream  develops  its  characteristic 
meanders,  vortex  shedding  and  ring  generation  capabilities  (Robinson, 
1971).  The  less  defined  and  broader  current  beyond  the  Grand  Banks  is 
then  known  as  the  North  Atlantic  Current.  Figure  1.1  is  a  schematic 
summarization  of  the  near  surface  currents  as  well  as  prominent 
topographic  features  as  presented  by  Watts  (1983). 

Typical  current  velocities  of  the  Gulf  Stream  between  Cape 
Hatteras  and  the  Grand  Banks  (hereafter  referred  to  as  the  Gulf 
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Figure  1.1     Schematic  Representation  of  Western  North 

Atlantic  Surface  Currents  (from  Watts,  1983) 


Stream)  are  on  the  order  of  1-2  knots.   The  dynamics  of  this  area  of 
the  Gulf  Stream  System  are  being  simulated  and  analyzed  in  this  study. 

B.   FORMATION  OF  EDDIES  IN  THE  GULF  STREAM 

The  presence  of  rings  associated  with-  the  Gulf  Stream  has  been 
documented  and  studied  as  early  as  1793  by  Jonathan  Williams,  grand 
nephew  of  Benjamin  Franklin  (Richardson,  1983).  Advancement  in 
techniques  of  their  location,  tracking,  and  documentation  as  well  as 
advances  in  the  physical  understanding  of  their  dynamics  have  led  to 
several  theories  concerning  their  formation.  These  theories  include 
dynamic  instabilities  (barotropic  and  baroclinic)  of  the  Stream  and 
topographic  forcing. 

The  rings  associated  with  the  Gulf  Stream  have  either  a  warm  or 
cold  core.  A  cold  core  ring  is  formed  by  Gulf  Stream  meanders  looping 
to  the  south  of  the  stream  as  depicted  in  Figure  1.2  (Richardson, 
1983,  as  adapted  from  Fuglister,  1972).  The  stems  of  the  meander 
merge,  trapping  a  core  of  slope  water  originally  located  to  the  north 
of  the  Stream  (Fuglister,  1972)  in  a  cyclonic  ring.  Similarly,  warm 
core  rings  form  anticyclones  to  the  north  of  the  jet  in  the  slope 
water  regions  and  have  as  a  core,  waters  entrapped  from  the  Sargasso 
Sea. 

Size  and  distribution  of  these  rings  vary  both  spatially  and 
temporally.  Figure  1.3  gives  a  synoptic  distribution  as  presented  by 
Richardson,  et  al.  ,  (1978).  Depicted  is  the  IS'C  isothermal  surface. 
Contours  are  based  on  expendable  bathythermograph,  CTD  hydrographic , 
and  satellite  infrared  data  from  March  16  until  July  9,  1975.   Notable 
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Figure   1.2 


Formation  of  a  Cold  Core  Ring 


Cold  core  ring  forcing  as  a  closed  segment  of  the  streamy 
circulating  around  a  mass  of  slope  water  (adapted  from 
Fuglister,  1972,  by  Richardson,  1983). 


Figure  1.3 


Gulf  Stream  Ring  Distribution 


Distribution  of  rings  from  March  16  to  July  9,  1975  as 
represented  by  the  IS'C  isotherm.   Primary  contour  intervals  are 
100  m  and  are  represented  by  the  darker  contour  lines  (from 
Richardson,  et  al. ,  1978). 


in  this  view  is  the  fact  that  though  the  radii  of  these  rings  vary, 
they  are  of  comparable  width  to  the  Gulf  Stream  itself.  It  is  on  the 
basis  of  these  observations  that  the  experimental  parameters  of  this 
study  are  based. 

C.  PROPAGATION  CHARACTERISTICS  OF  RINGS  WITHIN  THE  GULF  STREAM  SYSTEM 
Theories  for  the  movement  of  rings  within  the  Gulf  Stream  System 
vary,  with  numerous  recent  studies  shedding  additional  light  on  these 
propagation  mechanisms.  Observationally ,  the  movement  of  rings  is 
known  to  vary  depending  on  ring  location  relative  to  the  Gulf  Stream 
as  well  as  to  prominent  topographic  features  such  as  the  New  England 
Seamounts  and  the  continental  slope  and  shelf.  While  topographic 
steering  may  be  of  importance  near  the  shelf  and  slope,  the  primary 
emphasis  of  the  numerical  analysis  of  this  study  is  instead  on  the 
ring  and  Gulf  Stream  interactions  away  from  topography. 

Rings  that  are  in  contact  with  the  stream  tend  to  move  downstream 
at  speeds  of  up  to  75  cm/s  (Richardson,  1980).  Rings  that  are  not  in 
contact  with  the  stream  tend  to  propagate  with  a  westward  tendency 
(varying  from  northwest  to  southwest)  with  typical  speeds  of  5  cm/s 
(Lai  and  Richardson,  1977),  Additionally,  cyclonic  rings  may  have 
multiple  interactions  with  the  Gulf  Stream  and  their  propagation  will 
evolve  from  eastward  when  attached  to  the  stream,  to  southward  as  the 
ring  is  detached  from  the  stream,  then  westward,  and  followed  by 
northward  and  eastward  as  the  ring  again  interacts  with  the  stream 
(Richardson,  1980).  Additional  observations  by  Fuglister  (1977) 
confirm  this  "clockwise"  propagation  scheme  for  cyclonic  eddies. 


Figure  1.4 


Gulf  Stream  Ring  Trajectories 


The  trajectories  of  all  rings  continuously  tracked  with  free- 
drifting  buoys  plus  one  (ring  4)  tracked  by  SOFAR  floats  (see 
Cheney  et  al . ,  1976).   The  mean  path  of  the  Gulf  Stream  is  shown 
by  shading  (from  Richardson,  1980). 


Figure  1.4  provides  results  from  various  cruises  and  studies  that 
illustrate  each  of  the  propagation  characteristics  provided  above. 

There  have  been  different  mechanisms  proposed  for  ring  movement. 
The  first  is  due  to  the  representation  of  a  vortex  as  the 
superposition  of  many  Rossby  waves  (Mied  and  Lindemann,  1979). 
Associated  with  this  beta  dispersion- induced  westward  propagation  is  a 
subsequent  meridional  component  of  propagation  proportional  to  the 
strength  of  the  nonlinear  terms  (Firing  and  Beardsley,  1976).  In  this 
nonlinear  case,  the  gradient  in  the  eastern  portion  of  the  eddy  is 
strengthened  through  differential  westward  propagation  of  Rossby  waves 
allowing  for  northward  (southward)  propagation  components  of  cyclones 
(anticyclones).  In  addition,  rings  are  thought  to  move  due  to 
advection  by  large  scale  mean  flows.  Westward  movement  is  attributed 
to  advection  by  the  5  cm/s  recirculation  regime  of  the  Gulf  Stream 
System  (Richardson,  1983),  while  eastward  movement  is  attributed  to 
advection  by  the  Gulf  Stream  itself  when  the  eddy  is  attached  to  with 
the  stream  (Richardson,  1980). 

An  additional  theory  has  been  recently  suggested  by  Stern  and 
Flierl  (1987).  They  represent  a  Gulf  Stream  ring  by  a  point  vortex 
and  the  jet  by  a  jump  discontinuity  in  vorticity.  Their  results, 
obtained  on  an  f  plane,  give  westward  propagation  for  anticyclones 
(cyclones)  on  the  north  (south)  side  of  the  jet.  The  interaction  of 
an  anticyclone  (cyclone)  with  the  cyclonic  (anticyclonic)  side  of  the 
jet  leads  to  mutual  advection  westward,  much  the  same  as  the 
interaction  of  point  vortices  of  opposite  sign.  The  argument  for  this 
tendency  is  linear  for  eddies  adjacent  to  (but  not  in  the  immediate 
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vicinity  of)  the  jet.  They  give  an  analytic  solution  for  eddy 
propagation  speed  which  is  inversely  proportional  to  eddy  and  jet 
separation  distance.  For  eddies  close  to  the  jet,  nonlinear  effects 
become  important  and  eddy  propagation  speeds  are  found  by  using  the 
method  of  quasigeostrophic  contour  dynamics. 

D.   PURPOSE  OF  THIS  STUDY 

The  primary  objectives  of  this  study  are  the  numerical  analysis 
and  determination  of  ring  propagation  mechanisms  as  associated  with  a 
high  velocity  oceanic  jet.  These  are  determined  from  numerical 
experiments  in  which  the  initial  spatial  separation  of  the  jet  and 
ring  in  relation  to  each  other  is  varied.  Results  are  obtained  using 
a  two  layer,  nonlinear,  primitive  equation  model.  Experimental 
analyses  are  made  using  both  cyclonic  and  anticyclonic  vortices  on 
both  the  north  and  south  sides  of  an  eastward  flowing  jet.  Both  beta 
and  f  plane  simulations  are  examined.  While  cyclones  to  the  north  of 
the  jet  and  anticyclones  to  the  south  are  obviously  not  naturally 
occurring  phenomena  in  the  general  sense,  they  are  nonetheless 
included  for  the  study  and  understanding  of  the  dynamics  of  the 
problem  at  hand.  Results  are  compared  to  observed,  analytical  and 
numerical  propagation  tendencies.  The  results  support  specifically  a 
theory  of  westward  propagation  tendencies  related  to  vorticity 
advection  due  to  the  nonlinear  effects  associated  with  the  shear  at 
the  jet  and  ring  interface,  as  well  as  propagation  due  to  azimuthal 
perturbations  of  the  ring  due  again  to  nonlinear  interaction  between 
the  jet  and  ring.  The  results  also  show  enhanced  propagation  due  to 
beta  and  nonlinear  effects  on  isolated  vortices. 

a 


Complete  understanding  of  the  propagation  mechanisms  and 
characteristics  of  the  mesoscale  vortices  associated  with  the  Gulf 
Stream  is  necessary  for  a  more  thorough  understanding  of  the  general 
circulation  of  the  ocean  in  its  entirety  and  for  the  real  time 
prediction  of  eddy  location  and  movement  for  operational  use  in  the 
USN.  The  latter  is  due  to  the  extreme  acoustic  properties  associated 
with  the  warm  and  cold  core  rings. 


10 


II.   NUMERICAL  TECHNIQUES  AND  MODEL  PARAMETERS 

A.   THE  NUMERICAL  MODEL 

1 .  Model  Equations 

Simulation  of  the  Gulf  Stream  is  accomplished  utilizing  a  two 
layer,  primitive  equation,  semi  -  implicit ,  numerical  scheme.  This 
model  has  been  used  in  numerous  ocean  and  mesoscale  circulation 
studies  including  Hurlburt  and  Thompson  (1980,  1982),  Smith  and 
O'Brien  (1983),  and  Smith  (1986).  The  semi-implicit  model  was 
initially  developed  by  Smith  and  O'Brien  (1983).  Linear  test  cases 
were  performed  as  in  Smith  and  Reid  (1982)  that  compared  favorably 
with  linear  analytic  solutions  for  verification  of  Rossby  dispersion 
characteristics  of  the  model. 

Motion  in  each  layer  of  the  model  is  governed  by  a  momentum 
equation: 

^  +  (v-Vi+V^.v)  Vi  +  kxfVi  =  -hiVPi+AhWi         (2.1) 
and  a  continuity  equation: 

^  +   V.Vi^-0  (2.2) 

dt 

with  i  representing  layer  index  (i-1  upper  and  i=2  lower).  All 
variables  and  notation  are  found  in  Table  1.  The  fluid  is  Boussinesq 
and  hydrostatic  and  the  density  in  each  layer  is  constant. 

2 .  Model  Domain 

A  rectangular  finite  difference  gridded  domain  is  used  in  this 
model  with  grid  resolution  (2Ax)  of  20  km.   The  grid  is  oriented  in 
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PARAMETER 


TABLE  1 

REFERENCE  STATE  MODEL  PARAMETERS 

SYMBOL 


VALUE 


E-W  Basin  Extent 

N-S  Basin  Extent 

Initial  Upper  Layer  Thickness 

Initial  Lower  Layer  Thickness 

Maximum  Basin  Depth 

Depth  of  Bottom 

Coriolis  Parameter  ((§40°  N  lat.) 

df/dy 

Gravitational  Acceleration 

Reduced  Gravitational  Ace. 

Horizontal  Eddy  Viscosity  Coef . 

Time  Step 

Jet  Maximum,  i  th  Layer 

Eddy  Maximum,  i  th  Layer 

Jet  Position 

Eddy  Position 

Gradient  Operator 

Laplacian  Operator 

Pressure  in  Upper  Layer 
Pressure  in  Lower  Layer 


h 

1100  km 

Lx 

800  km 

hi 

1000  m 

h2(x,y) 

4000  m 

(h, 

+  ^z^max 

5000  m 

D(x 

y) 

5000  m 

fo 

0.94  X  lO--*  s"^ 

^0 

2.0  X  10-^^  m'^s- 

•  1 

g 

9.8  m  s-2 

_  f 

P2-Pi 

g 

'     P. 

Ah 

100  m^  s-^ 

At 

4200  s 

^ij 

100  cms"^ 

Vie 

100  cms"^ 

h 

variable 

Le 

variable 

V 

d       a 

ax   dy 

V2 

3x2  +  ay2 

Pi 

g(hi  +  hj) 

P2 

Px  -  &'\ 
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the  east-west  direction.  The  domain  encompasses  1100  km  in  the  x 
(east-west)  direction  and  800  km  in  the  y  (north-south)  direction. 
The  upper  layer  of  the  model  extends  to  1000  meters  while  the  lower 
extends  from  this  depth  to  5000  meters.  No  bottom  topography  is 
included . 

3 .   Boundary  Conditions 

Using  the  orientation  above,  free  slip  boundaries  are  used  on 
both  the  northern  and  southern  boundaries  with  no  flow  normal  to  the 
boundaries.  The  western  (upstream)  boundary  consists  of  a  prescribed 
inflow  velocity,  equal  in  both  the  upper  and  lower  layers.  The 
eastern  (downstream)  boundary  utilizes  a  radiation  condition 
(Carmerlengo  and  O'Brien,  1980),  in  which  flow  is  advected  out  of  the 
basin  at  speed  Ax/At  when  detected  adjacent  to  the  boundary.  In 
addition,  a  frictional  absorption  layer  is  implemented  on  this 
boundary  in  which  the  Laplacian  lateral  friction  coefficient 
(Aj^)  in  the  governing  equations  is  increased  from  250  m^s'^  (utilized 
elsewhere  throughout  the  domain)  to  1000  m^s'^  in  the  final  50  km  of 
the  domain  along  its  eastern  boundary. 

This  final  condition  is  obviously  not  representative  of  the 
conditions  within  the  Gulf  Stream  itself,  but  is  nonetheless  utilized 
in  order  to  prevent  back  radiation  of  small  scale  vorticity 
disturbances  into  the  active  areas  of  the  model  domain.  Various  model 
simulations  showed  that  the  east/west  dimension  of  the  domain  is 
extensive  in  comparison  and  the  dynamics  of  the  jet/ring  interactions 
were  not  affected  by  this  absorption  band. 
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B.   NUMERICAL  EXPERIMENTS 
1 .   Reference  State 

In  order  to  better  facilitate  the  isolation  of  ring 
propagation  mechanisms,  a  variety  of  model  simulations  and 
simplifications  were  made.  The  first  major  simplification  concerns 
the  stability  of  the  jet  itself.  The  Gulf  Stream  is  both 
barotropically  and  baroclinically  unstable.  The  mechanics  of  its 
meanders  and  the  continuous  variability  of  its  jet  axis  therefore 
impose  a  continuously  varying  distance  between  the  jet  and  any 
existing  ring  in  its  vicinity.  This  in  turn  imposes  additional 
variables  into  the  propagation  mechanisms  associated  with  jet/ring 
interactions  making  cross  correlation  of  variables  and  therefore 
isolation  of  specific  interactions  and  propagation  mechanisms  a 
somewhat  formidable  task.  Since  the  intent  of  this  study  is  the 
isolation  and  identification  of  these  propagation  mechanisms  at  and 
near  the  jet  interface,  a  barotropically  stable  jet  was  used 
throughout  this  study.  This  removes  from  the  problem  many  potentially 
complex  variables  due  to  jet  meandering.  It  also  simplifies  the 
identification  of  propagation  mechanisms  and  interface  processes  and 
provides  for  the  study  of  these  mechanisms  based  primarily  on  jet/ring 
proximity.  Accomplishing  the  identification  of  the  mechanisms  in  this 
manner  allows  for  the  establishment  of  a  reference  data  base  for 
future  use  in  more  realistic  numerical  simulations. 

The  next  simplification  made  in  this  study  is  the  absence  of 
the  southwestward  recirculation  currents.  Various  sources  and  studies 
list  westward  ring  advection  by  this  flow  as  a  valid  propagation 
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mechanism  since  the  magnitude  of  the  mean  return  circulation  is 
similar  to  that  of  typical  ring  propagation.  Worthington  (1976) 
documents  recirculation  speeds  of  5  cm/s  which  are  similar  to  the 
propagation  speeds  of  both  cyclonic  and  anticyclonic  rings.  In  this 
study,  all  simulations  are  conducted  without  any  southwestward  mean 
flow  with  the  intent  of  identifying  other  mechanisms  that  could 
possibly  duplicate  observed  propagation  speeds  and  directions. 

An  additional  simplification  involves  the  use  of  an  f  plane. 
Experiments  are  duplicated  on  both  f  and  beta  planes  for  the  further 
isolation  of  the  propagation  mechanisms  of  the  rings.  Whereas  beta 
effects  are  known  to  produce  westward  propagation  of  a  ring  (Warren, 
1967  and  Flierl,  1977),  an  intent  of  this  study  is  to  determine 
whether  other  propagation  mechanisms  are  also  valid,  in  particular  the 
point  vortex  theory  as  proposed  by  Stern  and  Flierl  (1987). 

Utilizing  the  aforementioned  simplifications  the  basic  state 
then  consists  of  a  two-layer,  geostrophically  balanced  Gaussian  jet  of 
the  form, 

h2(y)  =  H2  +  Ai[l-exp(-yV2Lj2,)]  (2.3) 

with  the  horizontal  scale  (Lp  equal  to  the  e-folding  width  scale  of 
the  jet.  Additionally,  a  two-layer,  geostrophically  balanced, 
Gaussian  vortex  is  included  as  represented  by, 

h2(x,y)  =  H2  +  A[exp(-RV2L%)]  (2.4) 
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h^(x.y)  -  H, 

with  R^  -  (x^  +  y^)  and  representing  the  radial  distribution  of  the 
ring. 

The  jet  e-folding  scale  is  set  at  45  km  while  the  eddy  e- 
folding  scale  is  set  at  40  km.  Figure  2.1  is  the  most  basic  case  of 
the  model  simulation  and  exemplifies  the  stability  of  the  jet  from  its 
initial  state  at  day  0  to  its  final  state  in  the  model  simulation  at 
day  36.  The  lines  within  the  figure  are  lines  of  constant  surface 
height  anomaly  and  are  representative  of  streamlines  of  the  barotropic 
current.  The  increments  on  the  x  and  y  axis  of  Figure  2.1  and  all 
subsequent  model  output  fields  represent  275  and  200  km  respectively. 
Unless  otherwise  specified  the  x  axis  is  1100  km  and  the  y  axis  is  800 
km  for  all  model  output  fields. 

The  initial  jet  and  ring  velocities  are  set  at  100  cm/s  in 
each  layer.  In  the  jet  the  100  cm/s  velocity  is  a  constant  inflow, 
while  that  of  the  ring  is  a  100  cm/s  gradient  balanced  initial 
condition  based  on  height  anomaly.  Figure  2.2  is  an  example  of  a 
model  field  output  of  velocity  contours  showing  the  comparable 
velocities  in  the  upper  and  lower  layers. 

For  the  reference  state  shown  previously,  and  for  all 
subsequent  model  output  field  plots  the  surface  height  is  contoured  in 
centimeters,  velocities  are  contoured  in  cm/s  and  potential  vorticity 
in  m"^  s'^ . 
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(a) 


(b) 


Figure  2.1 


Reference  Jet  Height  Ratio 


Height  fields  at  (a)  day  0  and  (b)  day  36.   Contours  are  labeled 
in  centimeters. 
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(a) 


Figure  2.2 


(b) 

Barocropic  VelociCy  Fields 


Velocity  fields  on  day  3  in  (a)  upper  layer  and  (b)  lower  layer 
Contours  are  labeled  in  cm/s . 
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2 .   Variation  of  Parameters 

As  the  focus  of  this  study  involves  vortex  propagation 
mechanisms  due  to  interaction  with  the  jet  interface,  the  major 
parameter  variation,  and  indeed  the  only  variable  used  other  than  f 
vs.  beta  planes,  is  the  vortex  and  jet  locations  in  the  initial  state 
of  each  model  simulation.  The  major  spatial  variations  utilized  are 
termed  "far"  and  "near"  fields.  In  far  field  cases,  the  vortex  is 
initially  located  at  distances  farther  than  3  jet  widths  from  the  jet. 
In  near  field  cases,  the  vortex  is  in  the  vicinity  of  the  jet  (i.e. 
outer  vorticity  and  height  anomaly  contours  of  the  jet  and  vortex 
intersect).   Any  other  location  is  termed  the  "mid"  field. 

The  experiments  were  divided  and  subdivided  as  follows.  The 
first  major  divisions  will  be  the  beta  plane  and  f  plane  simulations. 
Within  each  of  these  categories  fall  experiments  in  which  the  relative 
vorticity  of  the  vortex  is  of  the  same  sense  as  the  vorticity  shear  of 
that  side  of  the  jet  (for  example  a  cyclone  located  on  the  north  side 
of  the  jet) .  The  length  of  the  simulations  varies  from  21  to  45  days 
with  the  vast  majority  extending  to  36  days.  A  brief  description  of 
each  experiment  follows. 

a.   Beta  Plane  Simulations 

(1)  Experiment  Numbers  BTli3  and  hT26  (Initial  Beta  Plane 
Experiments) .  In  these  experiments  an  anticyclone  and  a  cyclone  with 
velocity  and  height  anomalies  the  same  as  in  all  other  jet  inclusive 
experiments  are  placed  on  a  beta  plane  without  the  presence  of  a  jet. 
The  purpose  of  this  was  to  determine  vortex  propagation  speeds  due 
solely  to  westward  Rossby  wave  propagation  and  nonlinear  meridional 
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effects.  This  excludes  the  effects  of  any  possible  propagation 
enhancement  or  retardation  due  to  the  presence  of  the  jet.  Initial 
and  final  states  are  given  in  Figure  2.3. 

(2)  Experiment  Numbers  BT35  and  BT4/3  (North  Side  Anti- 
cyclones) .  An  eastward  flowing  barotropic  jet  is  introduced  with  an 
anticyclone  to  the  north  of  the  jet  for  the  purpose  of  determining  jet 
influences  on  vortex  propagation.  Experiment  number  BT3^  is  the  far 
field  simulation  and  experiment  number  BT4^  a  mid  field  case.  Figure 
2.4  shows  the  height  anomaly  contours  and  exhibits  the  initial 
conditions  of  these  two  cases.  Figure  2.4  is  also  indicative  of  the 
initial  eddy/jet  separation  of  all  other  far  and  mid  field 
experiments . 

(3)  Experiment  Numbers  BT55.  and  BT6/3  (North  Side 
Cyclones) .  The  signs  of  the  rings  are  reversed  to  cyclones  and  cases 
of  far  and  near  fields  simulated.  This  near  field  case  was  simulated 
in  order  to  determine  if  the  mutual  effects  of  the  jet  and  ring 
interactions  could  actually  overcome  the  northwest  propagation 
tendencies  due  to  nonlinear  and  beta  effects  and  achieve  downstream 
advection  as  suggested  by  Stern  and  Flierl  (1987). 

(4)  Experiment  Numbers  BT7i9  and  BTSg  (South  Side 
Cyclones) .  These  experiments  are  as  essentially  mirror  images  of 
experiment  numbers  BT3)9  and  BT4^  with  cyclonic  rings  of  an  opposite 
sense  placed  to  the  south  of  the  jet  with  the  same  eddy  jet  separation 
distance.  Experiment  number  BT8/9  is  a  far  field  experiment  with 
experiment  number  BT7/9  being  a  near  field  case. 
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(a) 


day  0 


day  36 


^    ^ 

\ 

/ 

^ 

(b) 


day  0 


day  36 


Figure  2.3 


Isolated  Eddy  Cases 


Upper  layer  height  fields  for  (a)  an  anticyclonic  at  day  0  and 
day  36  and  (b)  a  cyclone  at  day  0  and  day  36. 
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(a) 


(b) 


Figure    2.4  Experimencs    BT3/3   and   BT4/3 


Height   fields    for    che    initial    state    of    (a)    BT3/3    (far    field)    and 
(b)    BT4/?    (mid   field) 
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(5)  Experiments  Numbers  BT95  and  BT10i9  (South  Side  Anti- 
cyclones) ■  The  rings  are  reversed  in  sign  from  BT7y9  and  BT8y9  and  are 
a  mirrored  duplicate  of  experiment  numbers  BT5/3  and  BT6/3 .  Experiment 
BT9^  is  a  mid  field  case  and  experiment  number  BTlO/9  is  an  extreme 
near  field  case. 

b.   f  Plane  Simulations 

The  f  plane  experiments  were  conducted  with  similar 
initial  conditions  as  that  of  the  beta  simulations.  These  experiments 
are  studied  and  discussed  in  greater  analytical  detail  than  the  beta 
experiments  and  are  therefore  listed  in  tabular  form  (Table  2)  for  a 
more  extensive  parameter  identification. 

Parameter  variables  listed  include  initial  location  L(0) 
of  the  eddy  relative  to  the  jet  edge.  Nonlinearity  (Q)  is  quantified 
as  a  maximum  velocity  ratio  (V^^j^)  to  the  maximum  Rossby  wave  speed, 
(^R^(j).  R(0)  is  a  non-dimensional  distance  between  the  jet  and  eddy 
and  is  included  for  direct  comparison  with  the  results  of  Stern  and 
Flierl  (1987).  Figure  2.5  is  a  schematic  of  representative  model 
parameters.  Also  in  this  figure  "C"  and  "AC"  represent  the  cyclonic 
and  anticyclonic  shear  associated  with  the  northern  and  southern  edges 
of  the  jet  respectively. 
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Figure  2.5     Schematic  of  the  Model  Domain  and  Initial  Eddy  and 
Jet  Configuration 
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III.   RESULTS  AND  ANALYSIS 

A.   BETA  PLANE  EXPERIMENTS 

1 .   Verification  of  Model  Output 

Various   beta   plane   experiments   were   conducted   and   are 

discussed  first  for  the  purpose  of  illustrating  ring  propagation 

associated  with  a  stable  jet  in  the  absence  of  a  recirculation  region. 

The  results  in  this  study  are  compared  to  previous  isolated  eddy 

numerical  studies  such  as  that  of  Mied  and  Lindemann  (1979).    In 

addition,  isolated  ring  propagation  is  compared  to  ring  propagation  in 

the  presence  of  a  jet  in  order  to  determine  if  in  fact  there  is  a  jet 

influence  on  the  propagation  tendency.   The  primary  purpose  here  is  to 

determine   which   "real   world"   results   may   be   associated   with 

barotropically  stable  beta  plane  simulations.   With  this  established, 

the  f  plane  simulations  can  then  be  conducted  in  order  to  isolate  in 

detail  what  additional  mechanisms  of  eddy  propagation  occur  due  solely 

to  the  jet/eddy  interactions.   By  using  this  approach,  time  dependent 

jet  spatial  variations  due  to  barotropic  and  baroclinic  instabilities, 

westward  beta  propagation  and  meridional  nonlinear  self  advection,  as 

well  as  recirculation  advection  are  systematically  removed  from  the 

equations  with  only  a  basic  state  remaining.    Thus  only  the  most 

rudimentary  aspects  of  jet/ring  interactions  remain  in  the  f  plane 

model  output. 

Mied  and  Lindemann  (1979)  conducted  numerical  experiments 

using  a  primitive  equation,  beta  plane  model  of  a  flat  bottom  two 
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layer  ocean  with  a  rigid  lid  imposed.  Their  experiments  were  of  upper 
ocean,  dispersing  and  pure  barotropic  eddies.  Their  results  indicated 
that  westward  propagation  is  associated  with  beta  effects,  and  with 
meridional  propagation  due  to  nonlinear  effects.  The  rate  of 
meridional  and  zonal  propagation  is  increased  with  increasing  current 
stre-gth  (i.e.  increased  nonlinear  self  advection) .  This  is 
attributed  to  the  increased  advection  of  planetary  vorticity  from  the 
north  (south)  of  the  eddy  to  either  side  of  the  cyclone  (anticyclone). 
As  the  strength  of  the  eddy  is  weakened  (numerically  through  viscous 
effects)  the  eddy  will  tend  to  turn  towards  the  west.  Their  results 
for  barotropic  eddies  indicated  that  meridional  propagation  varies 
from  2.5  to  9  km/day,  with  zonal  rates  in  the  2  to  3  km/day  range. 
The  results  of  experiment  BTl^  (Table  3)  are  not  inconsistent  with 
these  values . 

2 .   Opposite  Sign  Cases 

Ring  trajectories  for  experiments  in  which  the  eddy  is  of 
opposite  sign  as  that  of  the  shear  vorticity  of  the  jet  edge  (as  in 
figure  2.5)  are  shown  in  Figures  3.1a  and  3.1b.  These  and  all  other 
trajectory  plots  represent  30  days.  In  3.1a  an  isolated  anticyclone 
is  included  (BTl^)  along  with  a  mid  and  far  field  case  in  which  a  jet 
is  included  in  the  domain  to  the  south  of  the  eddy.  Translational 
velocity  averages  (in  9  day  increments)  for  each  experiment  are  given 
in  Table  3.  In  these  and  all  other  propagation  plots,  the  jet  axis 
for  vortices  to  the  north  of  the  jet  is  285  km.  For  the  cases  of 
vortices  to  the  south  of  the  jet,  the  jet  axis  is  at  555  km. 
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TABLE  3 
TIME  AVERAGED  BETA  PLANE  VELOCITIES  (NORTH  SIDE  ANTICYCLONES) 


EXPERIMENT 


DAY 

ZONAL  VELOCITY 

MERIDIONAL  VELOCITY 

(KM/DAY) 

(KM/DAY) 

9 

-2.0 

-3.8 

18 

-4.5 

-6.0 

11 

-1.0 

-5.2 

9 

-2.4 

-4.0 

18 

-1.8 

-8.3 

11 

4.2 

-1.6 

9 

-2.0 

-5.1 

18 

.5 

-6.0 

11 

5.2 

3.0 

BT1;9 


BT3^ 


BT4^ 
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Figure    3 . 1 


^  Plane  Vortex  Trajectories  (Opposite  Sign  Cases) 


30  day  trajectories  for  (a)  anticyclones  to  the  north  of  the  jet 
with  an  isolated  anticyclone  reference  (BTl^)  and  (b)  cyclones  to 
the  south  of  the  jet  with  an  isolated  cyclone  reference  (BT2^) 
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It  is  readily  apparent  from  both  Table  3  and  Figure  3.1a  that 
the  jet  has  little  influence  initially  on  the  ring  in  both  the  mid  and 
far  field  cases.  The  initial  propagations  are  southwest  as  in  the  no 
jet  case.  As  the  ring  approaches  the  jet,  propagation  direction  is 
reversed  as  the  ring  interacts  with  the  jet.  A  similar  phenomenon 
holds  true  for  the  zonal  propagation  tendencies  with  westward 
velocities  being  retarded  as  the  ring  approaches  the  jet  and 
velocities  being  reversed  to  eastward  as  the  ring  comes  into  contact 
and  interacts  with  the  jet.  The  results  here  are  mirrored 
duplications  of  the  clockwise  propagations  of  cyclones  to  the  south  of 
the  jet  as  observed  and  described  by  Richardson  (1980).  The  eddy 
motion  in  figure  3.1b  does  in  fact  resemble  these  observations  to  a 
large  extent  (see,  for  example,  ring  numbers  2,  3,  8  and  9  in  Figure 
1.4.  The  fact  that  Figures  3.1a  and  3.1b  are  mirrored  duplicates  of 
each  other  exemplifies  the  symmetry  of  the  dynamics  of  the  eddy  jet 
interaction. 

3 .   Same  Sign  Cases 

Figure  3.2  shows  beta  plane  trajectories  for  simulations  in 
which  the  sign  of  the  ring  velocity  is  the  same  as  that  of  the  shear 
vorticity  on  that  side  of  the  jet.  Though  these  simulations  are  of  a 
more  hypothetical  nature  than  the  opposite  sign  cases  discussed  above, 
they  are  nonetheless  included  for  comparison  again  with  suggestions  by 
Stern  and  Flierl  (1987).  They  show  that  like  sign  point  vortices  are 
initially  attracted  to  the  jet  while  "winding"  the  jet  interface 
around  the  vortex.  The  mutual  influence  tends  to  advect  the  ring 
downstream  at  speeds  similar  to  that  of  the  shear  flow.   The  near 
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field  cases  in  Figure  3.2  exhibit  tendencies  that  agree  qualitatively 
with  their  results.  BT5y9  and  BT9)8  on  the  other  hand  are  initially 
displaced  further  from  the  jet  and  the  nonlinear  propagation 
tendencies  overcome  the  influence  of  the  jet  and  allow  the  rings  to 
propagate  at  directions  and  speeds  similar  to  that  of  the  isolated 
ring  cases  discussed  earlier.  The  eddy  propagation  associated  with 
the  eddy  jet  interaction  (excluding  beta  effects)  is  more  clearly 
illustrated  in  f  plane  experiments. 

B.   f  PLANE  EXPERIMENTS 
1 .   Same  Sign  Cases 

The  same  sign  cases  (BT71F,  BT70F,  BT65F  and  BT66F)  have 
trajectories  as  shown  in  Figure  3.3.  Without  the  countering 
influences  of  beta  and  nonlinear  self  advection  the  eddies  are  able  to 
propagate  downstream  at  rates  as  high  as  30  km/day.  These  values  are 
similar  to  those  determined  by  Stern  and  Flierl.  Minor  differences 
are  due  to  the  finite  radius  of  the  eddy  and  width  of  the  jet  whose 
maximum  speed  is  at  its  central  axis.  In  the  numerical  simulations 
the,  ring  is  unable  to  completely  interact  with  the  jet  inner  core  as 
in  the  point  vortex  experiments  by  Stern  and  Flierl.  Also  of  note  is 
the  initial  attraction  of  the  rings  to  the  jet.  This  again  is  in 
agreement  with  Stern  and  Flierl.  Experiments  BT66F  and  BT71F  (far 
field  cases)  did  not  run  long  enough  to  advect  downstream  as  did  the 
two  near  field  cases,  but  of  note  is  their  eastward  (downstream)  turn 
at  the  latter  stages  of  each  run. 
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2 .   Opposite  Sign  Cases 

The  opposite  vorticity  sign  f  plane  cases  are  the  crux  of  this 
study.  Beta  induced  westward  motion  (Mied  and  Lindemann,  1979)  is 
thus  eliminated  as  a  propagation  mechanisms.  Meridional  motion  due 
to  nonlinear  effects,  subsequent  to  Rossby  dispersion  induced 
distortions,  is  also  removed.  The  southside  f  plane  trajectories  (as 
depicted  in  Figure  3.4)  as  well  as  all  cases  listed  in  Table  2, 
exhibit  a  propagation  tendency  and  therefore  a  propagation  mechanism 
that  is  not  associated  with  the  aforementioned  mechanisms.  Eddy 
propagation  paths  for  these  cases  can  be  categorized  as  meridional  and 
zonal.  In  the  near  field  cases  (small  L(0)),  the  initial  path  is  a 
meridional  ejection  away  from  the  jet  interface,  followed  by  a 
westward  drift.  This  response  is  particularly  evident  in  the  south 
side  cases  BT21F,  BT38fF  and  BT34F.  The  motion  for  larger  L(0)  is  a 
zonal  motion,  opposite  to  the  jet  direction  (in  this  case  westward). 

Figure  3.5  shows  the  evolution  of  the  potential  vorticity  in 
the  upper  layer  for  small  L(0)  experiment  BT21F.  This  experiment 
exhibits  the  evolution  of  vorticity  and  is  representative  of  all 
other  experiments  in  this  respect.  As  shown  in  Figure  3.4,  this  ring 
propagates  initially  perpendicular  to  the  jet  axis  at  6.5  km/day. 
This  southward  motion  is  visibly  less  prominent  as  L(0)  is  increased 
in  other  experiments.  The  mechanism  for  this  meridional  motion  is 
similar  in  nature  to  that  of  Stern  and  Flierl:  close  examination  of 
the  vorticity  plots  of  Figure  3.5  depicts  an  entrapment  of  jet 
vorticity.  The  region  between  the  southern  edge  of  the  jet  and  the 
eddy  is  one  of  strong  anticyclonic  shear.    Anticyclonic  relative 
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eddy  is  one  of  strong  anticyclonic  shear.  Anticyclonic  relative 
vorticity  is  then  advected  by  the  eddy  away  from  the  jet  and  is 
represented  in  this  figure  by  a  patch  of  anomalous  anticyclonic 
vorticity  to  the  west  of  the  cyclone.  The  pairing  of  the  anticyclonic 
and  cyclonic  vorticity  results  in  a  mutual  advection/steering  process 
that  results  in  the  initial  southward  propagation.  A  comparison  of 
BT21F  and  BT38F  also  indicates  that  the  pairing  of  eddy  and  jet 
vorticity  occurs  for  stronger  eddies  and  jets.  BT21F  has  initial  flow 
velocities  of  50  cm/s  in  contrast  to  the  remaining  experiments  in  this 
s  tudy . 

The  jet/eddy  interface  and  presence  of  this  anomalous 
vorticity  patch  also  results  in  a  perturbation  in  the  vortex  azimuthal 
structure  which  further  enhances  this  southward  propagation.  Figure 
3.6  shows  the  velocity  field  of  experiment  BT34F  on  day  6.  A  strong 
azimuthal  mode  1  perturbation  exists  with  intensified  gradients  on  the 
western  side  of  the  eddy.  As  the  eddy  weakens  and  propagates  away 
from  the  jet  axis,  the  anticyclonic  vorticity  patch  is  recaptured  by 
the  jet  and  is  advected  downstream  by  the  jet,  as  in  the  like  sign 
cases  discussed  previously.  As  it  propagates  and  rotates  to  the  north 
of  the  cyclone,  the  axis  of  the  azimuthal  mode  1  also  rotates  and  the 
advection  of  the  ring  is  thus  turned  to  the  west.  Figure  3.7  shows 
the  surface  height  field  of  near  field  experiment  BT38F  and  exhibits 
this  advection  and  rotation  of  the  axis  from  day  6  to  day  24. 

Stern  and  Flierl  attribute  these  tendencies  to  a  "capture"  of 
the  jet  interface  by  the  eddy  (Figure  3.8).  Their  results  indicate 
the  occurrence  of  this  capture  when  the  separation  distance  of  the 
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Figure  3.6 


Velocity  Field,  Experiment  BT34F  on  Day  6 
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"Capture"  of  the  Current  Interface 


"Capture"  of  the  current  by  an  anticyclonic  vortex  initially  at 
R(0)  -  0.5  and  located  at  the  indicated  point  R(6)  -  1.2  at  time 
t  -  6.   L(x,  6)  is  partially  indicated  by  the  solid  circles,  and 
the  remaining  three  curves  are  for  t  -  5.4,  and  2,  respectively. 
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vorticity  is  then  advected  by  the  eddy  away  from  the  jet  and  is 
represented  in  this  figure  by  a  patch  of  anomalous  anticyclonic 
vorticity  to  the  west  of  the  cyclone.  The  pairing  of  the  anticyclonic 
and  cyclonic  vorticity  results  in  a  mutual  advection/steering  process 
that  results  in  the  initial  southward  propagation.  A  comparison  of 
BT21F  and  BT38F  also  indicates  that  the  pairing  of  eddy  and  jet 
vorticity  occurs  for  stronger  eddies  and  jets.  BT21F  has  initial  flow 
velocities  of  50  cm/s  in  contrast  to  the  remaining  experiments  in  this 
s  tudy . 

The  jet/eddy  interface  and  presence  of  this  anomalous 
vorticity  patch  also  results  in  a  perturbation  in  the  vortex  azimuthal 
structure  which  further  enhances  this  southward  propagation.  Figure 
3.6  shows  the  velocity  field  of  experiment  BT34F  on  day  6.  A  strong 
azimuthal  mode  1  perturbation  exists  with  intensified  gradients  on  the 
western  side  of  the  eddy.  As  the  eddy  weakens  and  propagates  away 
from  the  jet  axis,  the  anticyclonic  vorticity  patch  is  recaptured  by 
the  jet  and  is  advected  downstream  by  the  jet,  as  in  the  like  sign 
cases  discussed  previously.  As  it  propagates  and  rotates  to  the  north 
of  the  cyclone,  the  axis  of  the  azimuthal  mode  1  also  rotates  and  the 
advection  of  the  ring  is  thus  turned  to  the  west.  Figure  3.7  shows 
the  surface  height  field  of  near  field  experiment  BT38F  and  exhibits 
this  advection  and  rotation  of  the  axis  from  day  6  to  day  24. 

Stern  and  Flierl  attribute  these  tendencies  to  a  "capture"  of 
the  jet  interface  by  the  eddy  (Figure  3.8).  Their  results  indicate 
the  occurrence  of  this  capture  when  the  separation  distance  of  the 
point  vortex  from  the  jet  is  less  than  a  specified  distance  T  (defined 
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as  the  square  root  of  vortex  circulation  divided  by  the  potential 
vorticity  of  the  shear  flow).  For  quantitative  comparison  the  jet 
edge  vorticity  in  this  numerical  study  was  estimated  as  0.5  *  10"*s"^. 
This  value  was  also  that  of  the  eddy,  due  to  its  similar  Gaussian 
structure.  For  the  eddy  radius  here  T  -  71  km.  Nondimensionalizing 
as  in  Stern  and  Flierl,  R(0)  -  L(0)/T.  Table  2  includes  values  for 
the  R(0)  used  within  this  study. 

For  barotropic  point  vortices,  Stern  and  Flierl  show  little 
meridional  motion  for  R(0)  >  2.4.  The  results  of  this  study  (Table  2) 
indicate  meridional  motion  for  R(0)  <  2.4  which  is  substantially  less 
for  values  greater  than  2.4.  These  results  are  therefore  in  agreement 
with  Stern  and  Flierl  with  respect  to  the  distance  within  which 
nonlinear  eddy/ jet  interactions  and  subsequent  propagation  occur. 

Perturbation  and  orientation  is  also  evident  in  mid  to  large 
R(0)  experiments  indicating  that  violent  collisions  and  interactions 
between  the  eddy  and  jet  are  not  necessary  for  this  meridional  motion 
to  occur.  Upper  layer  potential  vorticity  for  R(0)  -  1.6  (Figure  3.9) 
illustrates  a  weaker  capture  of  vorticity  from  the  jet  than  is  evident 
in  BT21F.  Also  evident  in  the  height  fields  is  the  persistent 
orientation  of  the  azimuthal  mode  1  structure.  Azimuthal  mode  1 
distortions  associated  with  Rossby  wave  dispersion  of  eddies  on  a 
beta  plane  have  been  seen  to  lead  to  meridional  vortex  motion  in 
nonlinear  eddies.  The  distortion  here  is  such  that  cyclone  (anti- 
cyclone) self -advection  to  the  south  (north)  is  opposite  to  that  of 
dispersion  induced  distortion. 
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Figure  3.9     Experiment  BT34F 

Contours  of  upper  layer  (a)  potential  vorticity  and  (b)  surface 
height  anomaly 
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In  addition  to  the  meridional  motions  described  above,  Stern 
and  Flierl  also  propose  a  theory  of  predominant  antiparallel  (zonal) 
propagation  associated  with  the  larger  R(0)  values.  They  suggest  a 
westward  motion  related  to  linear  vorticity  dynamics  in  which  the  eddy 
induces  a  perturbation  on  the  jet.  The  perturbation  vorticity  is 
opposite  in  sign  to  that  of  the  vortex  vorticity.  The  resulting 
vorticity  anomaly  interacts  with  the  eddy,  producing  westward 
propagation.  The  larger  separation  distances  prevent  the  more 
extensive  nonlinear  "capture"  of  the  interface  as  previously 
described. 

Experiment  BT28F  illustrates  this  nonlinear  capture.  A  small 
perturbation  is  seen  on  the  outer  jet  contours  in  Figure  3.10.  The 
perturbation  eventually  advects  downstream  but  imparts  a  westward 
propagation  to  the  eddy.  In  this  experiment  the  eddy  averages  1.4 
km/day  westward  motion.  Stern  and  Flierl  obtain  a  dimensional 
velocity  of  12  km/day  for  the  same  initial  separation  conditions. 
Although  this  is  a  factor  of  10  times  greater  than  that  determined  in 
BT28F  and  other  similar  experiments,  experiment  BT28F  is  more  in 
keeping  with  observational  studies.  Again  this  is  attributed  to  the 
more  realistic  model  representation  of  the  Gulf  Stream.  It  is  also 
interesting  to  note  that  eddy  azimuthal  distortions  seen  in  small  L(0) 
do  not  occur  in  large  L(0)  experiments. 
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Figure  3.10    Height  Anomaly  Field,  Experiment  BT28F 
Contours  of  height  anomaly  at  (a)  day  6,  (b)  day  12 
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IV.   DISCUSSION  AND  CONCLUSIONS 

A.   JET  INDUCED  EDDY  PROPAGATION  TENDENCIES 

Azimuthal  perturbations  to  Gulf  Stream  rings  have  been  observed  in 
repeated  hydrographic  surveys  of  a  ring  as  well  as  in  satellite  data. 
Hydrographic  surveys  indicate  that  a  ring  can  undergo  oscillations 
between  axisymmetric  mode  0  and  mode  1  or  mode  2  perturbations. 
Spence  and  Legeckis  (1981)  determined  that  a  cyclonic  ring  was 
elliptical  while  interacting  with  the  Gulf  Stream,  with  its  major- 
minor  axis  rotating  with  the  mean  flow  of  the  vortex  (cyclonically) . 
The  ring  returned  to  circular  (mode  0)  when  further  from  the  stream. 
Numerous  other  studies  have  also  shown  this  tendency  to  return  to 
circular  upon  ring  separation  from  the  Stream. 

This  study  indicates  that  a  strong  azimuthal  mode  1  perturbation 
is  induced  on  an  isolated  eddy  by  interaction  with  the  jet  when  the 
eddy  is  initially  close  to  the  jet.  The  orientation  is  such  that  the 
eddy  is  intensified  on  the  western  side  if  R(0)  is  sufficiently  small 
enough  to  "capture"  shear  vorticity  from  the  same  side  of  the  jet.  In 
this  case  the  eddy  vorticity  pairs  with  this  anomalous  patch  of 
opposite  sign  vorticity.  The  subsequent  meridional  motion  is  then  due 
both  to  the  western  intensification  of  the  ring  and  also  to  the  mutual 
advection  of  the  vortex  pair. 

As  the  vortex  (in  this  case  a  cyclone  to  the  south  of  the  jet) 
moves  away  from  the  jet,  the  anticyclonic  vorticity  anomaly  is 
advected  zonally  along  the  jet  axis.   This  in  turn  rotates  the  axis  of 
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the  mode  1  vortex  anticyclonically  resulting  in  a  north  side 
intensification  and  westward  motion.  As  vortex  and  jet  separate 
completely,  the  vortex  evolves  through  mode  2  to  axisymmetric  in  a 
time  scale  consistent  with  the  results  of  McCalpin  (1987).  Coincident 
with  the  evolution  to  axisymmetry  is  a  decrease  of  the  meridional  and 
zonal  propagation  rates. 

In  larger  R(0)  experiments  (>2.5),  the  initial  interaction  between 
jet  and  eddy  is  weaker  and  the  mode  1  perturbation  is  absent.  The 
eddy  remains  axisymmetric  and  interacts  weakly  with  a  perturbation 
induced  on  the  edge  of  the  jet.  This  results  in  predominantly  zonal 
motion  opposite  to  the  stream  direction. 

B.   EDDY  COALESCENCE  WITH  A  JET 

Although  the  focus  of  this  paper  is  not  on  the  coalescence  of  an 
eddy  with  a  jet,  certain  conditions  appear  necessary  for  coalescence 
to  occur.  The  experiments  here  show  no  tendency  for  eddies  to 
coalesce  with  the  jet  when  the  jet  is  in  a  stable  zonal  configuration. 
Instead,  an  eddy  interaction  with  the  jet  allows  the  eddy  to  capture 
opposite  sign  shear  vorticity  associated  with  the  jet  edge,  causing 
vortex  pairing  and  eddy  ejection  from  the   jet. 

Eddy  interactions  with  an  unstable  jet  can  be  rather  different. 
Additional  experimentation  showed  that  eddies  can  coalesce  with  an 
unstable  jet  when  a  collision  occurs  between  an  eddy  and  a  meander 
with  vorticity  of  the  same  sign.  Figure  4.1  shows  experiment  RG99 
which  illustrates  a  merger  of  a  cyclone  with  a  southward  extending 
cyclonic  meander.  Coalescence  occurs  due  to  the  comparable  vorticity 
of  eddy  and  meander.   An  experiment  with  an  anticyclone  in  the  same 
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Figure  4.1 


Experiment  RG99^ 


Contours  of  (a)  upper  layer  potential  vorticity  (b)  lower  layer 
potential  vorticity  and  (c)  height  field 
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location  did  not  indicate  a  merger,  but  instead  the  anticyclone 
advected  rapidly  eastward  with  the  jet. 

C.  RECOMMENDATIONS  FOR  FUTURE  STUDIES 

This  modeling  work  provides  a  base  on  which  to  build  and  further 
develop  studies  in  this  area  that  would  better  simulate  the  Gulf 
Stream  eddy  propagation.  A  first  step  would  be  the  utilization  of  a 
baroclinically  stable  and  then  an  unstable  jet  to  further  advance  the 
solution  of  the  actual  jet  induced  eddy  propagation  tendencies.  Next, 
imposing  periodic  inflow  and  outflow  boundaries  would  enable  for  a 
reduction  in  the  Aj^  lateral  friction  term  and  thus  allow  for  longer 
experimental  runs.  A  final  step  would  be  the  implementation  of  real 
hydrographic  and  satellite  data  in  the  initialization  of  the  model. 
If  approached  systematically  as  suggested  here,  this  would  eventually 
provide  for  real  time  eddy  propagation  and  jet  meander  forecasting  for 
operational  usage  in  the  USN. 

D.  CONCLUSIONS 

The  results  of  this  study  indicate  that  an  isolated  vortex 
interacting  with  a  comparable  strength  stable  jet  can  acquire  a 
meridional  propagation  tendency  away  from  the  jet.  The  translational 
speed  associated  with  eddy/jet  interactions  is  dependent  on  the 
initial  eddy/jet  separation  distance  and  the  strength  of  the  eddy  and 
jet.  For  Gulf  Stream  velocities,  the  eddy  translation  speed  away 
from  the  jet  increases  with  decreasing  eddy/jet  separation  if  the 
initial  separation  distance  is  less  than  about  4  times  the  internal 
Rossby  radius.   For  eddies  initially  within  this  distance  from  the 
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jet,  the  meridional  speed  is  proportional  to  the  strength  of  the  eddy 
and  jet.  For  eddies  initially  outside  this  range,  westward 
propagation  can  occur  which  is  associated  with  weaker  eddy/jet 
interactions.  The  meridional  and  zonal  propagation  tendencies  are 
related  to  vortex  pairing  between  the  eddy  and  opposite  sign  vortlcity 
of  the  jet  edge,  as  is  found  by  Stern  and  Flierl  (1987)  where  point 
vortices  and  a  jump  discontinuity  jet  were  used.  In  strong 
interaction  simulations,  the  motion  is  also  associated  with  strong 
azimuthal  mode  1  perturbations  acquired  by  the  eddy  in  the  eddy/jet 
interaction.  These  experiments  suggest  that  a  strong  ring  cannot 
coalesce  with  a  jet  when  the  jet  is  in  a  stable  configuration.  In 
experiments  where  the  jet  is  dynamically  unstable  however,  coalescence 
is  possible  only  when  an  eddy  collides  with  a  meander  of  like  sign 
vorticity.  In  this  case,  the  above  propagation  tendencies  are  no 
longer  observed. 
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